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ABSTRACT 

Microscopic calculations of neutron matter based on nuclear interactions de- 
rived from chiral effective field theory, combined with the recent observation of 
a 1.97 ± 0.04 Mq neutron star, constrain the equation of state of neutron-rich 
matter at sub- and supranuclear densities. We discuss in detail the allowed equa- 
tions of state and the impact of our results on the structure of neutron stars, the 
crust-core transition density, and the nuclear symmetry energy. In particular, 
we show that the predicted range for neutron star radii is robust. For use in as- 
trophysical simulations, we provide detailed numerical tables for a representative 
set of equations of state consistent with these constraints. 

Subject headings: dense matter, equation of state, neutron stars 

1. Introduction 

Neutron stars, apart from being systems for investigating such diverse topics as theories 
of gravity and the interstellar medium, are unique laboratories for studying matter at high 
densities. In neutron stars, matter ranges from nuclei embedded in a sea of electrons at 
low densities in the outer crust, to increasingly neutron-rich structures in the inner crust, 
to the extremely neutron-rich uniform matter in the outer core, and possibly exotic states 
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of matter at high densities in the inner core (Haensel et al. 2006 Lattimer 2012). The 
theoretical understanding of nuclear matter and atomic nuclei over such a range of densities 
and isospin asymmetry is a current frontier in nuclear science. The equation of state (EOS) of 
dense matter is also a key ingredient in modeling neutron star and black hole formation and, 



in particular, the gravitational wave signal from mergers of binary neutron stars (Andersson 



et al. 2011 Bauswein & Janka 2012 Bauswein et al. 2012) and neutron star-black hole 



mergers (Lackey et al. 2012) is sensitive to it. Consequently, future searches with advanced 



LIGO and LISA are expected to provide information about the high-density EOS. 

In the past, the EOSs of dense matter most commonly used in astrophysical simulations 
have been based on phenomenological nuclear interactions [for a review see, for example, 



Heiselberg & Pandharipande (2000)]. Two-nucleon (NN) interactions are usually constructed 
to fit NN scattering data at low energies. In addition, three-nucleon (3N) interactions were 
introduced because, first, they exist on theoretical grounds and, second, NN interactions 
alone cannot reproduce the properties of nuclei and nuclear matter. The many-body problem 
based on those interactions is very challenging due to strongly repulsive forces at small 
relative distances, which lead to highly correlated wave functions and require nonperturbative 
many-body methods. 

In recent years, the development of chiral effective field theory (EFT), following the 



pioneering work of Weinberg (1990, 1991), has provided the framework for a systematic 
expansion for nuclear forces at low momenta, where nucleons interact by pion exchanges 
and short-range contact interactions whose parameters can be fixed on the basis of two- and 



few-body observables (Epelbaum et al. 2009). Chiral EFT explains the hierarchy of two 



three-, and weaker higher-body forces and provides estimates of the theoretical uncertainties. 



In a recent Letter (Hebeler et al. 2010), we have shown that microscopic calculations based 
on chiral EFT interactions constrain the properties of neutron-rich matter up to nuclear 
saturation density to a high degree. On the basis of laboratory experiments and theory, 
our knowledge of the EOS at densities greater than 1 — 2 times the saturation density is 
limited. However, information may be obtained from measurements of neutron star masses. 
In particular, the recent discovery of a neutron star with a precisely determined mass of 1.97± 
0.04 M Q ( Demorest et al.|2010 ), the heaviest to date, is extremely important, as it rules out a 
large number of EOSs based on exotic degrees of freedom like hyperons or deconfined quarks. 
Such constituents soften the EOS at high densities and, without fine-tuning, are generally 
incompatible with a neutron star of such a large mass. Further information may be obtained, 



for example, from modeling X-ray bursts and quiescent low-mass X-ray binaries ( 


Ozel et al. 


2010 


Steiner et al. 


2010 


Ozel et al. 


2012 


Steiner et al. 


2013 


Giiver et al. 


2013 


), but the 



neutron star properties deduced are more model-dependent than the direct mass constraint 
from a very heavy neutron star. 
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By extending our microscopic results for the EOS at low densities in a general way to 
higher densities, we showed in Hebeler et al. (2010) that it is possible to derive systematic 
constraints on the EOS and on the radii of neutron stars. The high-density extensions we 
used were only constrained by causality and by the heaviest observed neutron star at that 
time which has a mass of 1.65 M Q . In this paper, we present details of these calculations 
and generalize and improve our approach in several ways: First, we require that the EOS is 
consistent with the observation of a 1.97 M Q neutron star (Demorest et al. 2010). To inves- 
tigate the sensitivity to the possible future discovery of neutron stars of higher mass, we also 
consider a second case where the EOS supports a neutron star of mass 2.4 M & . Second, we 
generalize the microscopic neutron matter calculations. Our previous results were based on 
renormalization-group-evolved NN interactions plus the leading 3N interactions from chiral 



EFT. The renormalization-group evolution improves the many-body convergence (Bogner et 



aL]|2010 ), but introduces uncertainties because the evolution was limited to NN forces. Here, 



we show that calculations based on unevolved chiral EFT interactions are in good agreement 
with the previous calculations (see Section [2|. Third, we improve the way beta equilibrium is 
incorporated (see Section [3]) and include explicitly the crust EOS below the crust-core transi- 
tion density, which we calculate in Section |4j Finally, we generalize the piecewise polytropic 
extensions of the EOS to higher densities by allowing more density regions and also refine the 
step size of the variation of the polytropic parameters (see Section [5]). We present our results 
for the nuclear EOS and the structure of neutron stars in Section |6l This shows that the 
constraints are robust and not significantly altered by the generalizations and improvements 
of the microscopic calculations and EOS extensions. For use in astrophysical simulations, 
we construct in Section [7] three representative EOSs consistent with the constraints from 
nuclear physics and observations and provide detailed tables in Appendix [Bj 



2. Neutron matter 



Our microscopic neutron matter calculations are based on chiral NN and 3N interactions. 
Neutron matter presents a unique system in chiral EFT because only the long-range two- 



pion-exchange parts of the leading 3N interactions contribute (Hebeler & Schwenk 2010). 



This is because three neutrons cannot interact via point-like S-wave interactions due to the 
Pauli principle. Moreover, the leading one-pion-exchange 3N interaction does not contribute 
in neutron systems because of the particular spin-momentum structure of this interaction. 

The leading chiral 3N forces have been shown to give important contributions to the 



nuclear EOS and to properties of nuclei (Hammer et al. 2013). In particular, saturation of 



symmetric nuclear matter is driven by 3N forces (Hebeler et al. 2011 ). While 3N contributions 
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to the neutron matter EOS are smaller in size than for nuclear matter, they are still significant 



and crucial for predictions of observables like the nuclear symmetry energy [see Hebeler & 



Schwenk (2010) and Section 3 . In this work, we include only the leading 3N forces, because 



it is presently possible to include only these beyond the Hartree-Fock level ( Tews et al.|2013 ). 



For calculations of the neutron matter EOS, the theoretical uncertainty is dominated 
by the uncertainties in the low-energy couplings c% and C3, which determine the two-pion- 
exchange three-body interactions between neutrons, not by the many-body approxima- 
tions (Hebeler & Schwenk 2010). These low-energy couplings relate 7rN, NN and 3N in- 



teractions. Their determination from 7rN scattering is within uncertainties consistent with 
the extraction from NN scattering, leading to C\ = —(0.7 — 1.4) GeV" 1 and c 3 = —(3.2 — 



5.7) GeV 


- 1 ( 


3ntem & Machleidt 


2003 




Epelbaum et al. 


2005 




Rentmeester et al. 


2003 


MeiBner 


2ooe 


). Beyond the leading 3N forces, part of the subleading 3N forces are sim- 



•(0.7 - 1.4) GeV -1 and c 3 = -(2.2 - 4.8) GeV" 1 for our calculations. 



In general, nuclear forces depend on an intrinsic resolution scale A. Consequently, the 
nuclear Hamiltonian can be written in the form 



if (A) = T + Vnn (A) + V 3N (A) + V 4 n(A) + . . . 



(1) 



where T denotes the kinetic energy, Vnn the NN interactions, V3N the 3N interactions, 
etc. The renormalization group provides a powerful tool to systematically change the scale 



A, while low-energy observables are preserved (Bogner et al. 2010). The evolution to low 
momentum scales improves the many-body convergence due to a decoupling of low and high 
momenta in the Hamiltonian (iB ogner et al.|[2~010 ). In general, the renormalization-group 
evolution changes all terms in Equation (JT]). The consistent evolution of 3N interactions in 



momentum space is a complex task and has been achieved only recently (Hebeler 2012). 
We start from the chiral N 3 LO NN potential with A = 500 MeV of |Entem fc Machleidt 



(2003) and use the renormalization group to evolve this NN potential to low-momentum 



Bogner et al. 


2007 


Hebeler et al. 


2007 



scales A = 1.8 - 2.8 fm 

teractions derived from different initial potentials are very similar (Bogner etalpOO^ . This 
universality can be attributed to common long-range pion physics and phase-shift equiv- 
alence. As a consequence, results of many-body calculations are rather insensitive to the 
particular choice of the initial NN interaction. Because the leading chiral 3N forces are of 
long-range character in neutron matter, they are expected to be to a good approximation 
invariant under the renormalization-group evolution for these cutoffs. Therefore, we use the 
leading chiral 3N forces determined by c\ and C3 also with low-momentum cutoffs. 



For low-momentum cutoffs, the EOS of neutron matter and nuclear matter can be 
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Fig. 1. — (Color online) Energy per particle of neutron matter as a function of density 
n. The blue band is based on chiral NN and 3N interactions with a renormalization-group 
evolution to improve the many-body convergence. The range of the band is mainly due to 
uncertainties in 3N forces ( Hebeler Schwenk||2010 ). The dashed red lines present the range 
without the renormalization-group evolution. 



calculated with theoretical uncertainties in a perturbative expansion in which the Hartree- 



Fock approximation is the first-order term (Hebeler & Schwenk 2010 Hebeler et al. 2011). 



Our neutron matter results also include second-order corrections [for calculational details 

For low momentum cutoffs 1.8 mi -1 < A < 2.8 fm -1 



sec 



Hebeler & Schwenk (2010) 



we 



have checked that contributions from third-order particle-particle diagrams give only small 
contributions to the energy per particle, about 25keV for A = 1.8 fm -1 and 300 keV for 
A = 2.8 fm -1 at the saturation density n = 0.16 fm -3 . In addition, we find that our 
second-order results are independent of the resolution scale to a very good approximation, 
and at the saturation density the maximal variation is about 1 MeV. We also find that the 
results are insensitive to the single-particle spectrum used. All these findings indicate that 
neutron matter is perturbative for low-momentum interactions, and show that the theoretical 



uncertainties of the many-body calculation are small (Hebeler & Schwenk 2010). 



Figure [T] shows our results for the energy per particle of neutron matter up to the 
saturation density (using a Hartree-Fock spectrum). The blue band is based on chiral NN 
and 3N interactions, with a renormalization-group evolution for NN interactions. The width 



of the band is due mainly to the uncertainties of c\ and C3 in 3N forces (Hebeler & Schwenk 
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Fig. 2. — (Color online) Pressure P of neutron matter as a function of density n. The blue 
band is as in Figure [T] based on chiral NN and 3N interactions. For comparison, the shaded 



red band shows the Quantum Monte Carlo results of Gandolfi et al. (2012) with 3N forces 
fitted to a symmetry energy of 32.0 — 35.1 MeV [see Table I of Gandolfi et al. (2012) . 



2010). For comparison, the dashed red lines present the range based on unevolved chiral 



NN interactions plus the same leading chiral 3N interactions. The remarkable agreement 
indicates that neutron matter is, to a good approximation, also perturbative for chiral NN 
interactions [for a detailed study, see Kriiger et al. (2013)]. We also explicitly checked the 
size of the contributions at third-order in the many-body expansion. As expected, they are 
larger for unevolved chiral NN interactions due to the stronger coupling between low and 
high momenta in the Hamiltonian. However, they are still significantly smaller than the 



second-order contributions (Kriiger et al. 2013) 



In Figure [2] we present the uncertainty band for the pressure of neutron matter based on 
chiral NN and 3N interactions (with a renormalization-group evolution for NN interactions). 



For comparison, we also show the Quantum Monte Carlo results of Gandolfi et al. (2012) 



(shaded red band) based on the phenomenological Argonne t> 18 NN potential plus 3N forces 
fitted to a symmetry energy of 32.0 



The 



35.1 MeV [see Table I of |Gandolfi etaL| ( [2012] ) 
agreement of the results is remarkable, given that the Hamiltonian and the many-body 
methods are completely different. 
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3. Asymmetric nuclear matter and beta equilibrium 

We extend the microscopic results for neutron matter to matter containing both neutrons 
and protons. To this end we use for the energy per particle e of asymmetric nuclear matter an 
expression that interpolates between the properties of symmetric nuclear matter and neutron 
matter. For e we take the kinetic energy plus an expression for the interaction energy that 
is quadratic in the neutron excess 1 — 2x: 

e(n,x) = -T \x 5 / 3 + (l-x) 5 / 3 ](2n) 2 / 3 
5 L J 

- [(2a - Aa L )x(l - x) + a L ] n + [(2r] - 4>q L )x(l - x) + r] L ] n 1 , (2) 

where n = n/n and x = n p /n denote the density in units of the saturation density and the 
proton fraction, respectively. T = (3n 2 n /2) 2 / 3 h 2 /(2m) = 36.84 MeV is the Fermi energy 
of symmetric nuclear matter at the saturation density. Equation ^ does not include the 
contributions from rest masses. The corresponding result for the pressure P = n 2 de/dn is 

P(n, x) = - n T \x 5/3 + (1 - x) 5/3 l (2nf /3 
5 L J 

- [(2a - Aa L ) x(l - x) + a L ] n 2 + 7 [(2 V - Ar] L ) x(l - x) + r} L ] n 7+1 . (3) 

The parameters a,r],aL and tjl can be determined from the saturation properties of sym- 
metric nuclear matter combined with the microscopic calculations for neutron matter of 
Section |2} For 7 = 4/3 and empirical saturation properties of symmetric nuclear matter, 

e (n = l,x = 1/2) = -B = -16MeV and P(n = l,x = 1/2) = , (4) 

this results in a = 5.87, 77 = 3.81, and a reasonable incompressibility parameter 

d 2 e(n, x 



K = 9 



dn 2 



= 236 MeV. (5) 

n=X,x=l/2 



The precise value of K can be adjusted by modifying the exponent 7 in Equations ^ and 
However, as shown in Table [TJ the predicted range for the symmetry energy and its density 
derivative depend very weakly on the particular choice of 7 and the resulting value of K, 
because the leading density dependence of the symmetry energy is linear. 

The parameters and t]l are extracted from the calculated bands for the neutron 
matter energy and pressure of Figures [T] and [2j Our results are based on the blue bands 
(with renormalization-group evolution), unless stated otherwise. We have first checked that 
the parametrizations ^ and (|3| provide excellent global fits for the energy and pressure up 
to a density n\ ~ 1.1 Uq. To determine qj, and t]l, we sample their values systematically 
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Fig. 3. — (Color online) Allowed range for aj, and i]l of the parametrizations ^ and ^ fit 
to the saturation point of symmetric nuclear matter and to the calculated neutron matter 
energy and pressure. 



7 


K [MeV] 


S v [MeV] 


L [MeV] 


1.2 


210 


29.7- 32.8 


32.4- 53.4 


4/3 


236 


29.7- 33.2 


32.5 - 57.0 


1.45 


260 


30.1 - 33.5 


33.6 - 56.7 



Table 1: Predicted range for the symmetry energy S v and for the L parameter, which de- 
termines the density dependence of the symmetry energy. Results are given for different 7 
values, which lead to different incompressibilities K, but, as shown, the predicted ranges for 
S v and L depend very weakly on 7. 

and require that the resulting energy and pressure be within the uncertainty bands shown 
in Figures [T] and [2] for densities from 0.45 no to 1.1 no- This leads to the allowed range for 
ol and rjL shown in Figure [3j with correlated limits c%l = 1.18 — 1.59 and t]l = 0.64 — 1.11. 

The proton fraction x for matter in beta equilibrium is determined by minimizing, for 
a given nucleon density, the total energy per particle, Equation (j2l), plus the contributions 
from electrons and from the rest mass of the nucleons. This amounts to the condition that 
^n + fn-nC 2 = fip + rripC 2 + u e , where \i n and \i v are the neutron and proton chemical potentials 
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n = n Q 


X 


fin [MeV] 


Hp [MeV] 


min 


0.040 


54.2 


-58.0 


max 


0.053 


51.9 


-71.5 


n = no/2 




min 


0.030 


34.6 


-46.1 


max 


0.033 


34.3 


-48.7 



Table 2: Proton fraction x and chemical potentials /i n and jj, p in beta equilibrium for the 
saturation density n and for no/2. The rows marked "min" and "max" give the range of 
the uncertainty band. 



without the rest mass contribution, or equivalently 



de(n, x) 
dx 



+ fjL e (n, x) - (m n - m p )c 







(6) 



For an ultrarelativistic, degenerate electron gas, the chemical potential is given by [i e (n, x) = 
hc(3Ti 2 xnon) 1 ^ 3 . The allowed ranges for and n^ imply ranges for the proton fraction 
and the neutron and proton chemical potentials in beta equilibrium, which are given for 
the saturation density no and for no/ 2 in Table [2j In the calculations we neglected the 
difference between the neutron and proton masses (1.3 MeV), which is small compared with 
fi e ~ 100 MeV. These ranges provide anchor points for other equations of state. 

The parametrizations ^ and ^ also make it possible to reliably extract the symmetry 
energy S v and its density derivative L, 



S „ 



1 d e(n, x 



8 dx 2 



n=l,x=l/2 



and L 



3d 3 



e(n, x 



8 dndx 2 



n=l,x=l/2 



(7) 



The region for ai and rji translates into an allowed region for the symmetry energy S v and 



the L parameter shown in Figure HI after Lattimer & Lim (2012). In addition, we give in 



Table [T] the predicted ranges for S v and L for different 7 values, corresponding to different 
incompressibilities K = 210 MeV, 236 MeV and 260 MeV. The predicted range for 7 = 4/3 
nearly spans the ranges for the other 7 values. This demonstrates that the extrapolation ^ 
is robust and that the theoretical uncertainty due to the choice of 7 is very weak and clearly 
much smaller than the empirical bands shown in Figure |4j 

In Figure [4], we compare the S v and L region predicted by our neutron matter results 
with values extracted from other data ( Lattimer &: Lim||2012 ). It is striking that the neutron 
matter results lead to the strongest constraints. These agree well with constraints obtained 



from energy-density functionals for nuclear masses (orange band) (Kortelainen et al. 2010) 
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24 26 28 30 32 34 36 
S v (MeV) 



Fig. 4. — (Color online) Constraints for the symmetry energy S v and the L parameter 
following Lattimer & Lim ( |2012 ). The blue region shows our neutron matter constraints, 
in comparison to bands based on different empirical extractions (for details see text). The 
white area gives the overlap region of the different empirical ranges. 



and from the 208 Pb dipole polarizability (yellow band) ( Tamii et al.||2011 ). In addition, there 



is good agreement with studies of the Sn neutron skin (light blue band) (Chen et al. 2010) 



of isotope diffusion in heavy ion collisions (HIC, green band) (Tsang et al. 2009), and of 



giant dipole resonances (GDR, red band) (Trippa et al. 12008 Lattimer & Lim 2012). More- 



over, there is very good agreement with an estimate obtained from modeling X-ray bursts 



and quiescent low-mass X-ray binaries (shaded region, labeled 'Astrophysics') (Steiner et al. 



2010). Remarkably, the constraints from these analyses have a common area of intersection 



(white area), which overlaps within uncertainties with the constraints from microscopic cal- 
culations of neutron matter. Based on the comparison in Figure [2j there is also very good 
agreement with the S v — L correlation band obtained from the Quantum Monte Carlo results 



of Gandolfi et al. (2012) [see Lattimer & Lim (2012); . 
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4. Crust-core boundary 



The transition between the neutron star crust and a uniform state in the core is a 
first-order phase transition. However, the density jump across the boundary is small; con- 
sequently, a good estimate for the density at which the transition from crustal matter to 
uniform matter in the core takes place may be obtained by determining the conditions under 
which, on lowering the density, matter becomes unstable to formation of a small density 



modulation (Baym et al. 1971a; Pethick et al. 1995). The density at which this occurs 



provides a lower bound on the density of the uniform phase at which the transition occurs. 

The electron screening length in the crust is large compared with typical nuclear separa- 
tions, so it is a good first approximation to regard the electron density as remaining uniform 
when a proton density wave is formed. For matter to be stable to formation of a long wave- 
length density fluctuation, the energy density must increase when the density modulations 
are imposed subject to the condition that the total number of neutrons and the total number 



of protons remain constant. Our discussion is adapted from Baym et al. (1971a) with minor 



variations. One condition for stability when Coulomb and contributions to the energy from 
density gradients are included is 



v 



+ 2(4vre 2 /3) 1/2 - 0k% T > 



The quantity 



dfj. p (dfx p /dn n ) 



dn r , 



(dfi n /dn n ) 



(9) 



is an effective proton-proton interaction when the Coulomb interaction is neglected. The 
first term is the contribution when the neutron density is constant and the second term 
is an induced interaction due to exchange of neutron density fluctuations. The effect of 
inhomogeneities in the density distribution is described by the quantity 0, given by 



(5 1{Qpp -f- ^QnpC QnnC 



(10) 



and a more detailed discussion is given in Appendix A Here Qij = Bij/n in the notation 

d\x v jdn n 



of Baym et al. (1971a) and 



^ dfi n /dn n ' 
The Thomas-Fermi wave number &tf is given by 

4 e 2 
k 2 --— k 2 

^TF — fe K e i 

71 he 



11) 



(12) 



with k e = (37r 2 nx) 1//3 . Generally, the second term in Equation (|sj) , which is due to the 
nonuniformity of the electron density, is small compared with the first but we do not drop it. 



30 31 32 33 34 30 31 32 33 34 

S„ (MeV) S v (MeV) 



Fig. 5. — (Color online) Contours of the lower bounds on the density of the liquid phase 
(solid lines) and the pressure (dashed lines) at the crust-core boundary in the S v — L plane. 
The density is measured in terms of the saturation density n and the pressure in MeVfm -3 . 
In the left panel, Coulomb and density gradient terms are neglected (Q = 0) and in the right 
panel Q = 75MeVfm 5 . The blue areas are the allowed region in Figure [4} 



If Equation (J8J) and the condition d\i n jdn n > are satisfied, matter is stable to small density 
modulations. With decreasing density, for realistic equations of state, the first stability 
condition to be violated is Equation (pi), and instability first sets in when this becomes an 



equality, which corresponds to Equation (9.18) of Baym et al. (1971a) 



Figure [5] shows contour plots of the densities and pressures at which instability sets in 
for uniform matter in beta equilibrium, with Coulomb and density gradient contributions 
for Q = 75MeVfm 5 (right panel), and without these contributions, Q = (left panel), see 
Appendix [Aj We find transition densities around n = 0.55 — 0.625, and inclusion of the 
Coulomb and density gradient contributions lowers these by about 15% to n = 0.475 — 
0.55. These transition densities are somewhat smaller than those predicted by the FPS and 



FPS21 interactions in Pethick et al. (1995), and consequently neutron star models with the 



interactions used in this paper will have lower crustal masses and lower crustal moments of 
inertia. 



5. General extension 



We investigate the structure of nonrotating neutron stars by solving the Tolman-Oppen- 
heimer-Volkov equations. Since the central densities of neutron stars can significantly exceed 
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Fig. 6. — (Color online) Pressure P of neutron star matter as a function of mass density p. 
The left part compares the pressure band predicted from the neutron matter results of 
Figure [2] and incorporating beta equilibrium following Section [3] to the BPS outer crust 



EOS (Baym et al. 1971b Negele & Vautherin 1973). The right part illustrates the general 



piecewise polytropic extension of the neutron star EOS to higher densities. We take p crust = 
Po/2, pi = 1.1 po, and vary the polytrope parameters over the ranges 1 ^ Ti ^ 4.5, 1.5 p ^ 
P12 ^ 8 po, ^ T 2 ^ 8, pi2 ^ P23 ^ 8.5 po, and 0.5 ^ r 2 3 ^ 8 (see text for details). 



the regime of our neutron matter calculations, we need to extend the EOS beyond this 
density. To this end, we employ a general polytopic extension, where the pressure of neutron 
star matter is piecewise given by -P(p) = np T , with mass density p = ran (Read et al. 



2009 Hebeler et al. 2010). This is illustrated in Figure [6j At low densities we use the BPS 
crust EOS ( |Baym et al.|1971b Negele & Vautherin 1973) up to the transition density p cms t- 
Figure [6] shows that for densities po/10 to po/2 the BPS crust EOS lies within the band 
predicted for the pressure of neutron star matter based on the neutron matter results and 
incorporating beta equilibrium following Section [3j [see also Hebeler et al. (2010)]. Therefore 



our results are insensitive to the particular choice of p crU st within this region. In the following, 
we use p CT ust = Po/2, based on our results for the crust-core boundary given by Figure [5j The 
pressure from p cr ust to pi = 1.1 po is given by the band predicted by chiral EFT interactions. 
Our results are insensitive to the precise value of pi in the vicinity of saturation density p , 
so we have taken a conservative value for which the uncertainty band of the microscopic 
neutron matter calculations is reasonable. 



For the extension beyond pi, we use three polytropes with exponents Ti,r2 and T3, 



-14- 



which make it possible to vary the soft- or stiffness of the EOS in the density regions 1: pi ^ 
p $C p 12 ; 2: pi2 ^ p ^ P23, and 3: p ^ p 23 , respectively (see Figure [6]). For densities just 
above pi, the EOS is still rather well constrained by the neutron matter calculations. The 
band predicted for the pressure of neutron star matter at pi corresponds to values of T in 
the range 2.25 — 2.5. Therefore, we take a restricted range for the first polytropic exponent 
1.0 ^ Ti ^ 4.5. We vary the value of all Tj in steps of 0.5. At intermediate densities we 
allow for the possibility of a phase transition and take a broad range ^ T 2 ^ 8. Finally, 
for densities beyond p 2 3 we allow for 0.5 ^ ^ 8. We exclude the value for this density 
region in order to avoid artefacts connected with a phase transition up to arbitrarily high 
density. For the densities between polytropes, pi 2 and p 2 3, we take 1.5 po ^ P12 < P23 < Pmax 
with a step size of po/2. We will show in the next section that the maximal density for our 
suite of EOSs of neutron stars is p max ~ 8.3 p . 

The general polytropic extension leads to a very large number of EOSs, which cover 
all possible pressures in the grey region in Figure [7j Furthermore we emphasize that this 
strategy is very general. It is based on a well defined uncertainty band at nuclear densities 
and does not rely on assumptions about the nature of the constituents of neutron star matter 
and their interactions at higher densities. The values of Tj and pij are limited by nuclear 
physics and observation. In the following, we will demonstrate that the recent observation 



of a 1.97 ± 0.04 M neutron star (Demorest et al. 2010) puts rather tight constraints on the 



EOS at high densities and the radii of neutron stars. We note that our results agree with 



those of our first study (Hebeler et al. 2010), which used only two polytropes. This shows 
that the general extension is robust and the conclusions would not change significantly if 
additional polytropes were introduced to characterize the pressure at high densities. 



6. Constraints on the nuclear equation of state and neutron stars 

The piecewise polytropic extension described in the previous section is used to generate 
a very large number of equations of state that cover the pressure-density plane at higher 
densities. We solve the Tolman-Oppenheimer-Volkov equations for each of these EOSs and 
retain only those that fulfill the following two constraints: 

1. ) the speed of sound remains smaller than the speed of light for all densities relevant in 

neutron stars: v* (p) = ^dP/dS c ^ c, where E is the energy density. 

2. ) the EOS supports a neutron star mass M = M, the mass of the heaviest neutron star 

observed. 
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Fig. 7. — (Color online) Pressure P of neutron star matter as a function of the mass density p. 
The blue band at lower densities represents the pressure predicted by the neutron matter 
results of Figure [2] with beta equilibrium incorporated as described in Section [3] The grey 
area shows the pressure accessed by the possible polytropic extensions. The lighter blue 
band is the envelope of the general polytropic extensions that are causal and support a 
neutron star of mass M = 1.97 M (left panel) and M = 2.4 M & (right panel). We also 
show the individual EOSs within these bands: The lines denote EOSs with central densities 
p c ^ 2.5 p (yellow), for 2.5 p < p c ^ 5 p (orange), and for p c > 5p (red). 



We consider each EOS up to densities at which the maximal neutron star mass is 
reached or the EOS becomes acausal, whichever density is smaller. In Figure [7] we present 
the individual EOSs that fulfill both constra ints for two cases: M = 1.97 M & (left panel), 
the mass of the heaviest known neutron star (Demorest et al. 2010), and M = 2.4 M & (right 
panel), an estimated mass of the black widow pulsar B1957+20 ( van Kerkwijk et al.| 2011 ). 
However, since the uncertainties of the latter determination are large, the M = 2.4 M 
constraint should be considered as a hypothetical mass, which is used here to probe the 
sensitivity of our results to the constraint from observations. 

The blue band at lower densities in Figure [7] represents the pressure predicted for matter 
in beta equilibrium as described in Sections [2] and [3] The lighter blue bands give the allowed 
EOS range, which is the envelope of the allowed polytropes at higher densities. The color of 
the individual lines indicates the maximal central density of the individual EOS (see Figure 
caption). Clearly, the pressure accessed by the possible polytropic extensions (the grey area 
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Fig. 8. — (Color online) Constraints for the pressure P of neutron star matter as a function 



of mass density p compared to EOSs commonly used to model neutron stars (Lattimer & 
Prakash 2001). The color coding of the bands is as in Figure [7j 



in Figure |7j) is substantially reduced by causality and by the constraint from neutron star 
mass measurements. The higher the mass of the heaviest neutron star observed, the stronger 
the EOS band is contrained. In addition, we find that the maximal densities in neutron stars 
are limited: stiff EOSs with large polytropic exponents have smaller maximal densities (see 
yellow lines), which are strongly constrained by causality. Softer EOSs tend to have larger 
central densities. For M = 1.97 M we find central densities as high as ~ 8.3 po, and, for 
M = 2.4 Mq the densities reach only « 5.8 po- 

In Figure[8]we compare the EOS bands of Figure [7]with a representative set of EOSs used 
in the literature. This set contains EOSs calculated within different theoretical approaches 
and based on different degrees of freedom. For details and notation we refer the reader 
to Lattimer & Prakash (2001 ). We find that a significant number of EOSs are not compatible 



with the lower density band based on chiral EFT interactions. In addition, at higher densities 
only very few EOSs, including the variational EOSs based on phenomenomenological nuclear 
potentials flAkmal et aL]|1998[ ), AP3 and AP4 in Fi gure [8j are located within the uncertainty 
bands over the entire density range. 



100 



1000 

£ [MeV fm" 3 ] 



Fig. 9. — (Color online) Constraints on the pressure P of neutron star matter as a function 
of energy density £. The color coding of the bands is as in Figure [7| The solid and dashed 
lines mark the allowed EOS range using the causal extension (see text) for M = 1.97 M 
(solid lines) and M = 2.4 M & (dashed lines). The dotted line gives the causal limit P = £. 

Figure [9] shows the uncertainty bands for the pressure as a function of energy density £. 
These are the natural variables to study to what extent the causality constraint is responsible 
for the apparent softening of the EOSs at high densities. For comparison we show the causal 
limit P = £, represented by the dotted line. Furthermore, we generated a causal extension for 
each individual EOS that reaches the causal limit at energy density ^it. For this we choose 
P(S) = £— £iimit.+-P(£iimit.) for £ > £ii m ;t, which ensures that the energy density, pressure, and 
speed of sound are continuous at all densities, with speed of sound v s {£) = c for £ > £u m it- 
In this way we can extend all EOSs to arbitrarily high densities and probe the role of the 
causality constraint on the results. In Figure |9j the solid lines mark the allowed EOS range 
using the causal extension for M = 1.97 M , the dashed lines the range for M = 2.4 M . 
The comparison to the blue bands shows that the causal extension changes only slightly the 
upper pressure limit and leads to somewhat higher maximal possible densities in a neutron 
star. We also observe that the causal extensions have a stronger impact on the lower limit of 
the pressure uncertainty band. For EOSs in this region the speed of sound reaches the speed 
of light already for small neutron star masses. By employing causal extensions, the maximal 
neutron star mass for such EOSs can increase significantly and consequently more EOSs 
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Fig. 10. — (Color online) Central densities p c of neutron stars with masses 1.4 M 
(green points), 1.97 M (red points), and 2.4 M (orange points) based on the uncertainty 
bands of Figure uj for the two cases M = 1.97 M & (left panel) and M = 2.4 M & (right panel). 





M = 


1.97 M Q 


M = 


2.4 M 




min 


max 


min 


max 


pjpo (1.4 M ) 


1.8 


4.4 


1.8 


2.7 


Pc /po (1.97 M ) 


2.0 


7.6 


2.0 


3.4 


p c /p (2.4 M Q ) 






2.2 


5.4 



Table 3: Minimal and maximal central densities p c (in units of the saturation density po) of 
the neutron stars shown by the points in Figure [lOj 



fulfill the mass constraint, which explains the significant difference between the blue bands 
and the regions within the black lines. Along the causal extension we have no information 
about the baryon density and therefore it is not possible to describe the equation of state 
there in terms of a (density-dependent) polytropic index d In P/d Inn. 



Figure 10 shows the central densities of neutron stars with masses 1.4M , 1.97M , 
and 2.4 M based on the uncertainty bands of Figure [7j The results for the minimal and 
maximal central densities are given in Table [3] Since stiff EOSs along the upper limit of 
the pressure uncertainty band are not sensitive to the mass constraint, the minimal central 
densities are identical for the M = 1.97 M and M = 2.4 M cases. 



8 10 12 14 16 8 10 12 14 16 

Radius [km] Radius [km] 



Fig. 11. — (Color online) Neutron star mass-radius constraints based on the uncertainty 
bands for the EOS of Figure Q for M = 1.97 M (left panel) and M = 2.4 M (right 
panel). The blue regions give the radius constraints based on the neutron matter results 
with renormalization-group-evolved interactions. The red dashed lines mark the band with- 
out renormalization-group evolution (see Figure [T]) and the thick dotted lines are for the 
allowed EOS range with causal extension (see Figure [9]). The solid line gives the causal 



limit (Lattimer 2012) 



The uncertainty bands for the EOS directly translate into constraints for the radii of 



neutron stars. In Figure 11, we present the radius constraints obtained from the EOS 
bands of Figure [7, for the two cases M = 1.97 M (left panel) and M = 2.4 M (right 
panel). The mass-radius relationships for the individual EOSs were obtained by solving 
the Tolman-Oppenheimer-Volkov equations and from these an envelope was constructed. 



The blue regions in Figure [TT] show the radius constraints based on the neutron matter 
results with renormalization-group-evolved interactions. In addition, we show results without 
renormalization-group evolution (see Figure [T]) and for the allowed EOS range with causal 
extension (see Figure [9]). The results with causal extension are also based on renormalization- 
group-evolved interactions. 

The predicted radius ranges are given in Table [4] for a 1.4 M neutron star and for the 
heaviest known neutron star, with M = 1.97 M . The maximal radius is very robust and 
essentially independent of the details of the neutron matter calculation, the use of causal 
extensions, and the mass constraints. This can be understood intuitively. Very stiff EOSs, 
which determine the maximal radius contraint, lead also to large neutron star masses. Hence, 
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unevolved 


causal 


evolved 


unevolved 


causal 
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(1.4M ) 


10.0 


9.7 


8.8 


11.6 


11.5 


10.4 


-Rmax 


(1.4 M ) 


13.7 


13.9 


13.7 


13.7 


13.9 


13.7 


-^min 


(1.97 M ) 


9.6 


9.3 


8.6 


12.0 


11.8 


10.9 




(1.97 M e ) 


14.2 


14.4 


14.4 


14.2 


14.4 


14.4 



Table 4: Radius constraints for a 1.4 M neutron star and for the heaviest known neutron 
star with M = 1.97 M based on the results of Figure 11 (see text for details on the column 



labels). The rows give the minimum and maximum radii. 



the constraints M = 1.97 M and M = 2.4 M are always fulfilled. Furthermore, the central 
densities of such neutron stars are typically rather low (see Figure 10). For the stiff est EOSs 
we find the central density to be p c ~ 1.8 po for an M = 1.4 M neutron star, which is 
typical for many observed neutron stars (see Table p|, a density region which is still rather 



well constrained by chiral EFT interactions (Hebeler & Schwenk 2010 Weise 2012). A 
1.4 Mq neutron star with a radius significantly larger than R = 13.9 km would therefore be 
incompatible with constraints from chiral EFT interactions. 

The minimal radius is more sensitive to details of the EOS at higher densities and there- 
fore less well constrained. The limits based on evolved and unevolved nuclear interactions 
agree well, except for very light neutron stars which are more sensitive to small differences 
in the low-density part of the EOS. In addition, the minimal radius strongly depends on 
the mass constraint. This implies that the lower limit of theoretical neutron star radii 
can be systematically improved by future observations of heavier neutron stars. Since the 
maximal mass of a neutron star based on soft EOSs can be increased by causal extensions 
(see the discussion of Figure [9]), we find a reduction of the minimal radius for EOSs by 
about 1 km with causal extensions, as shown by the difference between the blue bands and 



the dotted lines. This is consistent with the lower radius limits of Lattimer (2012), giving 
fi min (1.4M ) = 8.1km and 9.1km for M — 1.97 M and M — 2.4 M , respectively. 



7. Representative equations of state 

The uncertainty bands for the EOS are an envelope of many individual EOSs. Therefore, 
the limits of the uncertainty bands are not given by one EOS in general (see Figure [7]). In 
order to explore the EOS bands in astrophysical applications, we present three representative 
EOSs that probe the central region and the extremes of the uncertainty band for M = 
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Fig. 12. — (Color online) Left panel: Pressure P of neutron star matter as a function of 
energy density £ for three representative EOSs in comparison with the uncertainty bands 
of Figure [7j Right panel: The corresponding neutron star mass-radius results for the three 



EOSs and the uncertainty bands of Figure 11 



1.97 Mr. 



©■ 



a soft EOS (green dashed line in Figure 12), which agrees well with the lower limit 



of the EOS band for low and medium densities. This EOS describes excellently the 
minimal radius over the entire mass range (see the right panel of Figure 12). At 



higher densities, the pressure then reaches values above the lower limit of the band, 
as a consequence of the mass constraint. The polytropic parameters of this EOS are 



Ti = 1.5, pvi = 2.5 p , T 2 = 6.0, p 23 = 4.0 p , T 3 = 3.0, and p E 



7.0 p . 



2.) an intermediate EOS (orange solid line in Figure 12), which is located in the central 
region of the EOS band and is also consistent with the M = 2.4 M & EOS band over 
the whole density range. The radius of neutron stars for this EOS is about 12 km for 
typical masses, and lies in the center of the predicted radius range. The polytropic 
parameters of this EOS are Ti = 4.0, p i2 = 3.0 p , T 2 = 3.0, p 23 = 4.5 p , T 3 = 2.5, 
and p max « 5.4 p . 



3.) 



a stiff EOS (red dash-dotted line in Figure 12), which follows closely the upper limit 
of the EOS band until the causal limit is reached. This EOS gives the largest neutron 
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star radii consistent with the constraints. The polytropic parameters of this EOS are 
T 1 = 4.5, P12 = 1.5 p , T 2 = 5.5, p 2 3 = 2.0 p , T 3 = 3.0, and p max = 3.3 p . 

Detailed values of the pressure, energy density, and energy per particle for these three EOSs 
as a function of number density and mass density are listed in Tables [5] and [6j with the BPS 
crust EOS used at low densities in Table [7} In addition, Tables [5] and [6] include values for 
the neutron star radius and mass at that central density. 



8. Conclusions and outlook 

Recent advances in nuclear theory combined with new astrophysical observations have 
systematically tightened the constraints on the EOS of neutron-rich matter over a wide 
range of densities. In this paper, we have presented constraints for the EOS and for the 
structure of neutron stars based on microscopic neutron matter calculations with chiral EFT 
interactions. At nuclear saturation density and below, the uncertainties are dominated by 
uncertainties in 3N forces. In the future, these uncertainties can be reduced by improved 
constraints on the low-energy couplings entering 3N forces and by systematic improvements 



in chiral EFT that incorporate higher-order many-body forces (Tews et al. 2013; Kriiger et 



al.|2013[ ) and A de grees of freedom explicitly. In addition, reduced uncertainties are expected 



when renormalizat ion-group-evolved 3N forces ( Hebeler|2012 Hebeler fc Furnstahl|2013 ) and 



nonperturbative many-body calculations are employed. 

The uncertainties of the neutron matter calculations directly translate into constraints 
for the symmetry energy S v and its density derivative L, which is related to the pressure of 
neutron matter. On the assumptions that the interaction energy depends quadratically on 
the neutron excess, this leads to S v = 29.7 — 33.5 MeV and L = 32.4 — 57.0 MeV. These 
ranges provides very tight constraints and are consistent within uncertainties with different 
empirical extractions of these parameters (see Figure [4]). In addition, the neutron matter 
calculations predict the crust-core transition density to be in the range (0.475 — 0.55) n 
when Coulomb and density gradient corrections are included. 

We extended the EOS to higher densities by employing a general piecewise polytropic 
extrapolation, as illustrated in Figure [6} This led to a very large number of EOSs, which 
cover a wide pressure range, including the possibility of a soft phase transition region. From 
the individual EOSs we selected those that remain causal at all densities relevant in neutron 
stars and are able to support a neutron star with mass 1.97 M , the heaviest known neutron 



star (Demorest et al. 2010). We also considered a second case where the EOS supports a 



neutron star with mass 2.4 M . Note that the constraints are in accord with our previous 
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work (Hebeler et al. 2010), where we used only two polytropes. Combined with the mi- 



croscopic neutron matter calculations, this provides tight constraints on the nuclear EOS 
at sub- and supranuclear densities, and rules out many model EOSs, as demonstrated by 
Figure |8j 

Based on the allowed EOS band, we predicted ranges for the radii and central densities 

as well as Tables [3] and [4]). For the constraint 

13.9 km (based on 



11 



of neutron stars (see Figures 10 and 

M = 1.97 Mq, the radius of a 1.4 M Q neutron star is predicted to be 9.7 
renormalization-group-evolved or unevolved chiral EFT interactions) with central densities 
up to 4.4 pq. Note that the maximum radius is determined by the causality constraint, so 
that the discovery of a heavier neutron star only affects the minimum radius. If a 2.4 M 
neutron star were discovered, this would imply a radius range 11.5 — 13.9 km for a 1.4 M 
neutron star, with central densities up to only 2.7 po- 

For use in astrophysical simulations, we have constructed three representative EOSs 
(soft, intermediate, and stiff), consistent with the constraints from nuclear physics and ob- 
servation. In addition to the nuclear physics improvements mentioned above and the mea- 
surement of heavier neutron stars, information on the radii of neutron stars will significantly 
tighten the EOS band at high densities. Currently, observational radius limits often have 
conflicting results. Studies of five photospheric radius expansion X-ray bursts by |Steiner 
et al. (2010, 2013) are in agreement with the present analysis, but other analyses (Ozel et 



al. 2010) of two of them (EX01745-248 and 4U1820-30) indicate radii that are too small 



to be compatible with our results. On the other hand, a study of the burster 4U1724-307 



(Suleimanov et al. 2011) yields a radius too large to be compatible. Using phase modeling 



of pulsar profiles, Bogdanov (2013) found a compatible 3cr lower limit to the radius of pulsar 



PSR J0437-4715 of 11.1 km, assuming the 1.56M lower limit for its mass found by Verbiest 



etaL] ( |2007j ). However, |Hambaryan et al.| ( |20TT] ) found R = (10 - 12) (M/M Q ) km for the 
pulsar RBS 1223 using the same technique, which is compatible only if M < 1.4M . Also, 



Leahy et al. (2011) found that the sources SAX J1808.4-3658 and XTE J1807-294 had radii 



which were compatible, but that XTE J1824-338 did not, being too large. A recent analysis 
( Guillot et al.|201~3~ ) of 5 quiescent neutron stars in low-mass X-ray binaries in globular clus- 
ters, in which it was assumed that the stars all had the same radius, determined the radius 
to be R = Q.llj'gkm to 90% confidence, which is marginally consistent with our results. 
Additional radius measurements from the LOFT ( Feroci et al.||2012 ) and NICER (Gendreau 



et al. 2012) X-ray missions are therefore eagerly anticipated. 
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A. Instability of the uniform phase 



Here we discuss details of the density-gradient contributions to the instability condition 
and also express the instability condition in terms of the variables n and x, as was done in 



Lattimer & Prakash (2007), rather than n n and n p . The contribution of density gradients 



to the energy density is expressed in the form (Baym et al. 1971a) 

£ g md = <5 PP (Vn p ) 2 + 2Q pn (Vn p ) ■ (Vn n ) + Qnn(Vn n ) 



(Al) 



where the coefficients Qij are in general functions of the neutron and proton densities, but 
will here be treated as constants. We follow Baym et al. (1971a) and take Q np = Q pn = 
2Qnn = 2Q PP , in which case the quantity (3 of Section [4] is given by 



/3 = 2Q„ n (l + 4C + C 2 



(A2) 



Indeed, realistic Skyrme density functionals such as Ska and SkM* satisfy this approximate 



relationship among the Qij values (Lattimer & Prakash 2007). It would be interesting to 
determine the Qij from modern energy density functionals (Erler et al. 2012a[b ). 



The quantity Q may be determined either from the surface energy of symmetric nuclear 
matter or from the surface thickness of symmetric nuclei. For symmetric nuclear matter, the 
energy density is given by 



£buik + £ g rad = n(e s (n) + mc 2 ) + Q(Vnf 



(A3) 



where n[e s (n) +mc 2 ] = £(n, x = 1/2), Q = 3Q nn /2 = ?>Q np /A and m = (m n + m p )/2. In the 
Thomas- Fermi approximation and for symmetric nuclear matter the optimum density profile 
for a plane surface is determined by minimizing the quantity f dz [n(e s (n) — fi ) + Q{n') 2 ] 
for a surface lying in the xy-plane (Bethe 1968). Here fiQ = —B is the chemical potential of 



symmetric nuclear matter at the saturation density and the prime denotes a derivative with 
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respect to z. This leads to the condition 



n[e a (n) - /i ] h T 

n=- ] j — ^ — ( ) 

when the nuclear matter is assumed to occupy the region with negative z. This equation 
defines f(n/n ), and, on integration, determines n(z). Note that f(n/n ) — > for both 
n — > and n — )■ n . The surface energy per unit area, the surface tension cr , is then given 
by 



a = 2I a ^QT n 3 , (A5) 

where 

I a = / « 0.1696, (A6) 

for the bulk energy per particle given by Equation (|2| with 7 = 4/3, a = 5.87, and 77 = 3.81 
[see details after Equation pi)]. The 90% — 10% surface thickness can be expressed as 



£90-10 

where 



°- 9 du 
0.1 /(«) 



^, (A7) 



4.776 . (A8) 



For Q = 75MeVfm 5 , one finds a w 1.15 MeVfm -2 and £90-10 ~ 2.7fm, both reasonable 
values although perhaps 5% too large for the surface thickness. Note that Hartree-Fock 
calculations of surface energies can differ by about 10% from the Thomas-Fermi approxi- 



mation (Ravenhall et al. 1972). However, we have no freedom to fit both observables with 



the adopted functional for the bulk energy per particle. With this value for Q, we find 
Q np = lOOMeVfm 5 and Q pp = Q nn = 50MeVfm 5 , which are the values used in the calcula- 
tions in Section HI 

We now comment briefly on the stability conditions when expressed in terms of the 



variables n and x (Lattimer & Prakash 2007), rather than n n and n p . To take care of 



the constraints on the numbers of neutrons and protons, it is convenient to work with the 
thermodynamic potential 

E = £- fi° n n n - fi° p n p , (A9) 

where and n p are the chemical potentials in the initial state. With this choice, the first 
order terms in an expansion of H in powers of deviations of the densities from those in the 



initial state vanish. With n n = n(l — x) and n p = nx, it follows from Equation (A9) that 



* 2s - 55? (fa » 2 + (H + "° " "°) SnSx + \w (fa)2 • (A10) 



-26- 



where the term with the chemical potentials comes from the nonlinear dependence of the 
neutron and proton densities on n and x. Since the baryon pressure is given by P = 
n 2 d(£/n)/dn and the difference of the neutron and proton chemical potentials by n [fi p — 
fi n ) = d£/dx, it follows that 

d 2 £ _ 1 dP d 2 £ 1 dP ,1 d 2 £ 

n dn 



1 dP 

n dx 



A*n + A* P , and 



n- 



d{fi p - fi n ) 



(All) 



dn 2 n dn dxdn n dx 2 dx 2 dx 

Therefore, since in the coefficients we may put fi n = fi n and /i° = fi p , the quadratic form, 
Equation (A10), may be rewritten as 

1 dP 



5 Z E 



2n dn 



1 dP 

n dx 



(5n) 2 + -^-5n5x+ " '"'"^ 



n d(fi, 
2 



dx 



(5x) 



(A12) 



The conditions for the quadratic form to be positive definite are that the diagonal terms be 
positive, 

dP „ , d(ji p -ji n ) 



> and 



> 



dn dx 
and that the determinant of the quadratic form be positive (Lattimer & Prakash 2007[ ) 

dPd(i2 p -fi n ) fdP^ 2 



(A13) 



dn 



dx 



- - 



\ dx 



> 



(A14) 



If the condition (A14) and one of the conditions (A13) are satisfied, the other condition (A13) 
is satisfied automatically. If Equation (A13) is satisfied, the inequality (A14) may be rewrit- 
ten as 



cP£_ 
dn 2 



dP 

dx 



dx 



> 0. 



(A15) 



d(fi p - fi n ) 

The first term may be regarded as a direct interaction between density fluctuations, without 
changes in the composition, while the second represents an induced interaction between 
density fluctuations due to changes in the composition. For the system to be stable, the 
total interaction consisting of direct and induced contributions must be positive. 
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Table 6. Continuation of Table [5} The maximal density for the stiff EOS is 

reached already in Table [5j 
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6.0 


902.0 


419.7 


1077 


183.0 


9.90 


1.92 




















6.1 


917.0 


441.1 


1102 


190.3 


9.88 


1.94 




















6.2 


932.1 


463.1 


1128 


197.8 


9.85 


1.95 




















6.3 


947.1 


485.8 


1154 


205.4 


9.83 


1.96 




















6.4 


962.1 


509.3 


1180 


213.1 


9.80 


1.97 




















6.5 


977.1 


533.6 


1206 


220.9 


9.77 


1.98 




















6.6 


992.2 


558.6 


1233 


228.9 


9.75 


1.99 




















6.7 


1007 


584.4 


1261 


237.0 


9.72 


1.99 




















6.8 


1022 


611.0 


1289 


245.2 


9.69 


2.00 




















6.9 


1037 


638.4 


1317 


253.5 


9.66 


2.01 




















7.0 


1052 


666.5 


1345 


261.9 


9.63 


2.01 





















intermediate 
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Table 7. 



Numerical data of the BPS E OS for the outer crust ( Baym et al. 1971b Negele & 
Vautherin 1973). The units are as in Table [§ 



n/n Q 




pc 2 




P 




£ 




2.956 10" 


14 


4.444 10" 


12 


6.303 10" 


2b 


4.385 10" 


12 


2.975 10- 


14 


4.472 10" 


12 


6.303 10- 


24 


4.407 10- 


12 


3.069 10- 


14 


4.613 10- 


12 


6.303 10- 


23 


4.547 10" 


12 


4.369 10- 


14 


6.568 10- 


12 


7.551 10- 


22 


6.472 10- 


12 


6.18710- 


14 


9.302 10- 


12 


8.737 10- 


21 


9.150 10- 


12 


1.700 10- 


13 


2.556 10" 


11 


1.061 10- 


19 


2.516 10" 


11 


7.937 10- 


13 


1.193 10" 


10 


3.632 10- 


18 


1.183 10- 


10 


4.331 10 _ 


12 


6.511 10- 


10 


1.186 10- 


16 


6.416 10" 


10 


3.934 10" 


11 


5.915 10- 


-9 


6.081 10- 


15 


5.825 10- 


-9 


9.881 10- 


11 


1.485 10" 


-8 


3.100 10- 


11 


1.463 10- 


-8 


2.483 10- 


10 


3.732 10" 


-8 


1.51710- 


13 


3.675 10" 


-8 


6.235 10- 


10 


9.373 10- 


-8 


7.183 10- 


13 


9.228 10" 


-8 


1.566 10" 


-9 


2.355 10" 


-7 


3.286 10- 


12 


2.319 10" 


-7 


3.934 10" 


-9 


5.91410" 


-7 


1.447 10- 


11 


5.825 10" 


-7 


9.881 10" 


-9 


1.485 10" 


-6 


6.088 10- 


11 


1.463 10" 


-6 


2.483 10" 


-8 


3.732 10" 


-6 


2.441 10- 


10 


3.676 10" 


-6 


3.125 10" 


-8 


4.698 10" 


-6 


3.282 10- 


10 


4.627 10" 


-6 


6.235 10" 


-8 


9.373 10" 


-6 


8.955 10" 


10 


9.233 10" 


-6 


1.244 10" 


-7 


1.870 10" 


-5 


2.392 10" 


-9 


1.842 10" 


-5 


2.482 10" 


-7 


3.732 10" 


-5 


6.278 10" 


-9 


3.676 10" 


-5 


4.953 10" 


-7 


7.445 10" 


-5 


1.625 10" 


-8 


7.337 10" 


-5 


9.881 10" 


-7 


1.485 10" 


-4 


4.166 10" 


-8 


1.464 10" 


-4 


1.244 10" 


-6 


1.870 10" 


-4 


5.453 10 


-8 


1.843 10" 


-4 


1.972 10" 


-6 


2.964 10" 


-4 


1.01710" 


-7 


2.922 10" 


-4 


3.125 10" 


-6 


4.698 10- 


-4 


1.890 10" 


-7 


4.631 10" 


-4 


3.934 10" 


-6 


5.91410" 


-4 


2.57710" 


-7 


5.830 10" 


-4 


4.952 10" 


a 
o 


7.445 10" 


-4 


3.143 10" 


7 


7.342 10" 


-4 


6.235 10" 


-6 


9.373 10" 


-4 


4.281 10" 


-7 


9.245 10" 


-4 


9.881 10" 


-6 


1.485 10" 


-3 


7.938 10" 


-7 


1.465 10" 


-3 


1.566 10" 


-5 


2.355 10" 


-3 


1.470 10" 


-6 


2.323 10 


-3 


2.482 10" 


-5 


3.732 10" 


-3 


2.722 10 


-6 


3.683 10" 


-3 


3.125 10" 


-5 


4.698 10- 


-3 


3.533 10" 


-6 


4.637 10" 


-3 


3.934 10" 


-5 


5.91410" 


-3 


4.80710" 


-6 


5.836 10" 


-3 


4.952 10" 


-5 


7.445 10- 


-3 


6.540 10" 


-6 


7.353 10" 


-3 


6.235 10" 


-5 


9.373 10" 


-3 


8.893 10" 


-6 


9.256 10" 


-3 


7.850 10" 


-5 


1.180 10" 


-2 


1.209 10" 


-5 


1.166 10" 


-2 


9.881 10" 


-5 


1.485 10" 


-2 


1.562 10" 


-5 


1.468 10" 


-2 


1.244 10" 


-4 


1.870 10" 


-2 


2.124 10" 


-5 


1.848 10" 


-2 


1.566 10" 


-4 


2.355 10- 


-2 


2.888 10" 


-5 


2.328 10" 


-2 



n/n 




pc 2 




P 




£ 




1.972 10" 


-4 


2.964 10" 


-2 


3.713 10" 


-b 


2.931 10" 


-2 


2.482 10" 


-4 


3.732 10- 


-2 


5.048 10" 


-5 


3.692 10- 


-2 


3.125 10- 


-4 


4.698 10- 


-2 


6.865 10" 


-5 


4.649 10" 


-2 


3.934 10" 


-4 


5.91410 - 


-2 


9.330 10- 


-5 


5.852 10" 


-2 


4.953 10" 


-4 


7.445 10- 


-2 


1.269 10- 


-4 


7.376 10" 


-2 


6.235 10- 


-4 


9.373 10- 


-2 


1.621 10" 


-4 


9.284 10- 


-2 


6.906 10" 


- 1 


1.038 10- 


-1 


1.805 10" 


-4 


1.029 10- 


-1 


7.850 10- 


-4 


1.180 10- 


-1 


2.053 10" 


-4 


1.169 10- 


-1 


9.881 10" 


-4 


1.485 10- 


-1 


2.791 10" 


-4 


1.473 10- 


-1 


1.244 10" 


-3 


1.870 10- 


-1 


3.630 10" 


-4 


1.855 10- 


-1 


1.608 10" 


-3 


2.41710- 


-1 


4.871 10- 


-4 


2.398 10- 


-1 


1.669 10- 


-3 


2.509 10- 


-1 


4.924 10" 


-4 


2.488 10- 


-1 


1.95410" 


-3 


2.937 10- 


-1 


5.212 10" 


-4 


2.91710 - 


-1 


2.469 10" 


-3 


3.712 10- 


-1 


5.678 10" 


-4 


3.688 10- 


-1 


2.974 10" 


-3 


4.471 10- 


-1 


6.135 10- 


-4 


4.443 10- 


-1 


3.633 10" 


-3 


5.461 10- 


-1 


6.759 10" 


-4 


5.427 10- 


-1 


4.464 10" 


-3 


6.711 10- 


-1 


7.601 10- 


-4 


6.673 10- 


-1 


5.491 10" 


-3 


8.255 10- 


-1 


8.731 10- 


-4 


8.207 10" 


-1 


6.250 10- 


-3 


9.396 10- 


-1 


5.944 10" 


-4 


9.390 10- 


-1 


2.500 10- 


-2 


3.758 




5.799 10- 


-3 


3.762 




5.000 10" 


-2 


7.517 




1.166 10" 


-2 


7.530 




7.500 10" 


-2 


11.27 




2.085 10- 


-2 


11.30 




1.000 10" 


-1 


15.03 




3.21610" 


-2 


15.08 




1.250 10- 


-1 


18.79 




4.515 10" 


-2 


18.86 




1.500 10" 


-1 


22.55 




5.961 10" 


-2 


22.64 




1.750 10- 


-1 


26.31 




7.544 10" 


-2 


26.42 




2.000 10" 


_ i 


30.07 




9.260 10" 


_2 


30.21 




2.250 10- 


-1 


33.82 




1.11010" 


-1 


34.00 




2.500 10- 


-1 


37.58 




1.308 10" 


-1 


37.79 




2.750 10- 


-1 


41.34 




1.51810" 


-1 


41.58 




3.000 10" 


-1 


45.10 




1.742 10" 


-1 


45.38 




3.250 10- 


-1 


48.86 




1.980 10" 


-1 


49.18 




3.500 10- 


-1 


52.62 




2.231 10" 


-1 


52.98 




3.750 10- 


-1 


56.37 




2.49710" 


-1 


56.78 




4.000 10" 


-1 


60.13 




2.77710" 


-1 


60.58 




4.250 10- 


-1 


63.89 




3.073 10" 


-1 


64.38 




4.500 10- 


-1 


67.65 




3.385 10" 


-1 


68.19 




4.750 10- 


-1 


71.41 




3.71310" 


-1 


72.00 




5.000 10" 


-1 


75.17 




4.054 10" 


-1 


75.81 
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